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It is proposed that critical balance — a scale-by-scale balance between the linear prop- 
agation and nonlinear interaction time scales — can be used as a universal scaling con- 
jecture for determining the spectra of strong turbulence in anisotropic wave systems. 
Magnetohydrodynamic (MHD), rotating and stratified turbulence are considered under 
this assumption and, in particular, a novel and experimentally testable energy cascade 
scenario is proposed for low-Rossby-number rotating turbulence. It is argued that in 
neutral fluids, the critically balanced anisotropic cascade provides a natural path from 
strong anisotropy at large scales to isotropic Kolmogorov turbulence at very small scales. 
It is also argued that the kj 2 spectra seen in recent numerical simulations of low-Rossby- 
number rotating turbulence are analogous to the k ± spectra of the numerical MHD 
turbulence in the sense that they could be explained by assuming that fluctuations are 
polarised (aligned) approximately as in an inertial wave (Alfven wave for MHD). 



1. Introduction 

Ability to support both linear waves and nonlinear interactions is ubiquitous in natural 
systems. Wave turbulence is, therefore, a very generic situation in such systems when 
dissipation coefficients are small and ener gy injected at some s ystem-specific scale has 



to be dissipated at much smaller scales (Za kharov et a/.l 119921 ). Theory of turbulence 



is concerned with the ways in which the energy is transferred from large (injection) to 
small (dissipation) scales and, consequently, with the structure of the fluctuations in the 
intervening scale range. 

A common property of many such systems is the presence of some mean field that 
introduces a special direction. Examples are plasmas embedded in a mean magnetic field, 
rotating fluids and stably stratified fluids with a mean temperature or density gradient 
(in real systems usually in the direction of gravity) . Both linear and nonlinear physics is 
affected by the mean field: turbulent fluctuations in such systems tend to display a high 
degree of anisotropy. The typical wave numbers parallel and perpendicular to the special 
direction associated with the mean field often satisfy k\\ <C k±, while the wave dispersion 
relation is of the form 

U) = fc||U(fc±). (1.1) 
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For Alfven waves in magnetohydrodynamics (MHD), v — va, the Alfven speedJH for 
inertial waves in rotating fluids, v = 2Q/kj_, where Q is the rotation frequency. Our 
arguments will apply directly to these two cases; in stratified turbulence, fcii ^> k±, so 
the roles of k\\ and k± are reversed and certain adjustments to the general argument will 
be needed — they are explained in section [5l 

The dispersion relation of the type (11. 1| implies that waves propagate primarily in the 
parallel direction: indeed, the parallel and perpendicular group velocities are vu = v(k±) 
and v± — (k\\/k±)k±v'(Jc±) <C V\\. If the nonlinearity is of the fluid type, u ■ Vit, where 
u is the fluid velocity, then ku <C implies that nonlinear interactions are primarily 
perpendicular: u ■ Vu ~ u± ■ Vj_u, so the nonlinear decorrelation time is given bjQ 



r^L ~ k±u±(k±), (1.2) 



where u±(k±) is the characteristic velocity fluctuation amplitude corresponding to the 
wave number k±. Note that incompressibility V • u = and fcii -C k± imply that the 
perpendicular motions are individually incompressible, Vj_ • u± = (see appendix [A|l . 

For anisotropic wave systems, Kolmogorov-style dimensional theory alone does not fix 
the scalings of the energy spectra. Indeed, assuming a local (in scale) energy cascade and 
hence a scale-independent energy flux e, 



k±E(k±) ~ ui(fcx) ~ er(k±), (1.3) 

where E(k±) is the one-dimensional perpendicular energy spectrum and r(k±) is the 
"cascade time" corresponding to the characteristic wave number k±. In the absence of 
waves or anisotropy, it is dimensionally inevitable that r(k) ~ tnl(&)> whence E(k±) ~ 
e 2 / 3 fc~ 5 / 3 , the Kolmogorov spectrum. However, with waves and anisotropy, two additional 
dimcnsionless ratios arise: ku/kx and wtnl ~ k^v/k±u±, which measure the strength 
of the anisotropy and the relative time scales of the linear propagation and nonlinear 
interaction (equivalently, the relative size of the fluid velocity and the wave phase speed) . 
The spectrum can, as far as dimensional theory is concerned, be an arbitrary function of 
these two ratios, both of which can have some nontrivial scaling with k±. 

Clearly, an additional physical assumption is necessary to fix the scalings. In strong 
MHD turbulence, it is known as the critical balance (CB) and states that the ch aracter- 
istic l inear and nonlinear times are approximately equal at all scales: u> ~ ( Higdonl 



Il984t iGoldreich fc Sridharlll995l ). We propose that CB be adopted more generally as a 
universal scaling conjecture for anisotropic wave turbulence. This removes the dimen- 
sional ambiguity in determining the cascade time, so we may set r ~ tnl and recover 
the Kolmogorov spectrum, 

E(k±) ~ e 2 / 3 fc; 5/3 . (1.4) 

The CB itself, u> ~ , then gives us a relationship between the parallel and perpendic- 
ular scales: 

~ [TNL^xM^r 1 ~ £ 1/3 [v^T 1 kT- (1.5) 
Let us first give a general argument in favor of this idea (section ^ and then discuss 
three examples: MHD (and, more generally, plasma) turbulence, from whence these ideas 
originate (section [3]), rotating turbulence, for which we propose a novel energy cascade 
scenario (section [4j, and stratified turbulence (section [5]). Note that in section l4~5l we 
discuss the extension to the rotating turbulence of the concept of polarisation alignment 



f More generally, low-frequency waves in magnetised plasma s satisfy the gyrokinetic disper- 
sion relation, whi ch is also of the form (|Howes et aZ.ll2006ft . 

J See section 1331 for a discussion of when this is valid and what happens when it is not. 
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(also originating from MHD turbulence; see lBoldvrevll2006h. which may help interpret 
thefcj 2 spectra reported in numerical simulations ( Mininni et al. 2009; Th iele fc Miillerl 
2009). The section on rotating turbulence is the main part of this paper, while the MHD 
and stratified cases are discussed only briefly to emphasise what appears to be universal 
nature of some of the scaling arguments involved. In our concluding remarks (section [6]), 
we will also mention a few other examples of CB emerging as a general physical principle 
in systems that are different from the anisotropic wave type discussed here. 



2. Why anisotropic turbulence is neither weak nor two-dimensional 

In turbulent wave systems, if the fluctuation amplitudes at the injection scale are so 
small that wtnl 3> 1, the nonlinearity can be treated perturbativel y and what is known a s 



weak turbulence theory emerges as a controlled approximation (jZakharov et all 119921 ). 



In anisotropic wave systems, it typically predicts a turbulent casca de primarily in k 



(because the nonlinearity is primarily perpendicular), at constant fen (jGaltier et aLll200 



)0; 



Galtier 20031 ). While the analytic calculations can be quite involved, the basic result can 



be recovered in a simple nonrigorous way. If wtnl 3> 1, nonlinear interactions between 
wave packets result in small perturbations of the amplitudes 6u± ~ (^Tnl) These 
perturbations can be assumed to accumulate as a random walk and then the cascade time 
r is by definition the time that it takes the cumulative perturbation to become comparable 
to the amplitude itself: n^'^Sux ~ u±, where n ~ ra is the number of interactions. This 
gives t ~ ^t"nl ancl , using equations (|1.HI1.3[) . we get the one-dimensional perpendicular 
energy spectrurrfj] 

E(kx)^(ek ll ) 1 / 2 [v(kx)] 1/2 kZ 2 (2.1) 
The weak turbulence condition under which equation (|2.1|) is valid is then wtnl ~ 
fc||u(/c±)e _1 ' /3 fc ± 2 ^ 3 ^> 1. Unless v(k±) increases sufficiently fast with k±, the nonlinearity 
becomes stronger with increasing k± compared to the linear propagation and the weak 
turbulence condition is broken at k± given by equation (|1.5|l . Thus, the weak turbulence 
cascade drives itself into a critically balanced state. 

The opposite limit is a pure two-dimensional (2D) state: fcii is assumed so small that 
wtnl *C 1 and wave propagation is neglected. As generically happens in 2D, the energy 
cascade should then be inverse, from larger to smaller k±. As k± decreases, tnl becomes 
longer, so the 2D approximation, wtnl *C 1, is eventually broken and CB is reached, at 
which point the turbulence is again three-dimensional (3D). Thus, both from the weak- 
turbulence limit (small amplitudes) and the 2D limit, the turbulence naturally evolves 
towards a state of CB. 

There exists another argument, which is independent of the assumption of inverse 
cascade and suggests that 2D motions are fundamentally unstable. Consider two perpen- 
dicular planes separated by some distance. The motions in each plane will decorrelate 
on the time scale tnl- In the parallel direction, information is transmitted by waves, 
so perfect correlation between the two planes required for a pure 2D state can only be 
sustained if a wave can propagate between them in a time shorter than tnl- Thus, for 
any given k±, there will be some parallel distance, fc^ 1 , given by the CB relation (|1.5|) . 
beyond which the motions will decorrelate and become 3D. Thus, an initially 2D pertur- 

f We emphasise that there is no parallel cascade and the energy injection rate e can be an 
arbitrary function of fen. Note that in some other anisotropic wave systems there can be also a 
cascade in the parallel direction. The weak turbulence sp ectra for su c h syst ems (historically the 
first example of anisotropic wave spectra) were found bv lKuznetsovl ll 19721 ). 
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bation will tend to a state of CE0. This process can be interpreted as an instability of 
the 2D motions with respect to Cherenkov-type emission of waves. 



3. MHD (Alfvenic) turbulence 

The ideas laid out above in a general form originate from considerations of MHD 
turbulence. In MHD, equation (| 1 . 1 1) describes Alfven waves with v — va = Bo/y/ATrp, 
where Bq is the mean magnetic fiel d and p t he density of the conducting fluid. Both 
weak-turbulence ( Galtier et aLlEoOoh and 2D ( Montgomery fe Turner 198ll ) theories for 
MHD turbulence have been proposed. By the general arguments given above, both will 
naturally evolve towards a CB state, k\\VA ~ kj_uj_, with a Kolmogorov spectrum (11. 4| 
and a scale-dependent anisotropy given by equation (|1.5[) : 



(3.1) 



Note that as k± increases, the turbulence becomes more anisotropic (ku/k± decreases). 

If the turbulence is weak at the injection scale, its spectrum is expected to be [see 
equation (|2.1[) ] 

E(k ± ) ~ ( ekwVA^kZ 2 (3.2) 

and there is no cascade in k\\ ( Galtier et aZ.ll200oh . From equation (|3.2j) . u TIOS ~ (ew^/fco) 1 / 4 
where ko is the wave number of energy injection (assumed isotropic), so e ~ M\v A ko, 
where Ma = u tu1s /va -C 1 is the Alfvenic Mach number. Then the wave number at 
which weak turbulence breaks down and the critically balanced cascade begins is, from 
equation (|3.ip . 

fc± c ~ e- 1/2 (fc||^) 3 / 2 ~ k Q M- 2 . (3.3) 

The anisotropy of Alfvenic turb ulence and even equation (13.11) appear to have been con- 
firmed by numerical simul ations (ICho fe Vis huiac bood: iMaron fc Goldreichll200lf ) and 
solar wind measurements ( Horburv et aZ.ll2008 ; Podestall2009l). The weak tur bulence spec- 
tral scaling (|3.2|) has also been checked numerically ( Perez fc Boldvrevll2008l ). The precise 
nature of the scaling in the CB regime remains somewhat mysterious: while the solar wind 

-5/3 



measurements support k. 



(e.g 



Bale et a/.N2005t iHorburv et aLll2008t iSahraoui et < 
20091). numerical simulations give spectra much closer to k, 3/2 |m aron fc Goldreicbl 2001: 
Mason et aLll2008l ) — a modified critical balance argument proposed bv lBoldvrevi (|2006l) 



to explain these results will be discussed in section [4751 where we will show how it can 
be adapted to the case of rotating turbulence. 

3.1. Plasma turbulence 

Beyond the MHD approximation, the gyrokinetic disp ersion relation for low- frequency 
waves in magnetised plasmas is also of the form (jl.lj) (jHowes et all 120061 ) . The general 
idea of a critically balanced cascade can be extended to various typ es of gyrokinetic 
turbulence, e.g., for plasma t urbulence below the ion Larm or scale (|Cho fc Lazarianl 
12004 ISchekochihin et a/.ll2009l) . Further details can be found in lSchekochihin et all (|2009h 
(Sec. 7); the important point to keep in mind is that generalising the argument proposed 
in the present paper to plasma systems requires correctly identifying the cascading quan- 
tity (not always kinetic energy) and the type of nonlinearity (not always u ■ Vu). 

In the case of the turbulence of kinetic Alfven waves (dispersive waves that replace the 
MHD Alfven waves below the ion Larmor scale), the scaling predictions resulting from 
the application of the CB conjecture to their dispersion relation (also of the form (|1.1| ) 



X This argument was suggested by S. C. Cowley. 
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have been confirmed numerically (jBiskamp et aLlll999UCho k Lazarianl l2004. 2009); the 
sub-Larmor-sc ale spectra measured i n the solar wind also appear to be consistent with the 
CB prediction ( Sahraoui et a/.ll2009l) . This is the first example of confirmed applicability 
of the CB conjecture beyond its original target of MHD turbulence. 



4. Rotating turbulence 

Are magnetic anisotropy and magnetised plasma waves a special case or can CB be 
adopted as a universal scaling conjecture? A key test of universality would be to find a 
critically balanced cascade in a purely hydrodynamic setting. We propose the following 
scenario for the rotating turbulence. 

4.1. Critically balanced rotating turbulence and restoration of isotropy 
The dispersion relation for inertial waves in a rotating incompressible fluid is uj — 2f2/c|| jk. 
Suppose that the energy is injected isotropically at some characteristic wave number ko 
and the Rossby number Ro = u rms fco/^ ^ lj i-e., the amplitudes at the injection scale 
are so low that turbulence is weak. Then the en ergy cascade will proceed anisotropically, 
at fc|| ~ ko and towards larger kj_ (lGaltierll2003h . When kj_ 3> ku , the dispersion relation 
takes the form (jl.ip with v(k±) = 2Q/k±; the weak turbulence spectrum is then given 
by equation (|2.ip [analogous to equation (|3.2[) ] : 

E(k±) ~ (efc tt) 1/2 fc; 5/2 . (4.1) 

Note that this implies w rms ~ (ett) 1 / 4 ^ 1 ^ 2 and so e ~ Ro 4 ^ 3 ^ 2 . As k± increases, the 
nonlinearity becomes stronger and CB is reached at a critical k± that can inferred from 
equation (|1.5|) by setting fc|| ~ ko [analogous to equation (|3.3p ]: 

k± c ~ £^ 1/5 (fc||r!) 3/5 ~ fc R O - 4/5 . (4.2) 

For k± > k± c , the turbulence is no longer weak, but it is still anisotropic and the 
cascade is critically balanced: the spectrum is given by equation (|1.4| . while equation 
(|1.5[) becomes 

fc,l ~ s^Q^k 5 / 3 . (4.3) 

This relation is qualitatively different from the MHD case [equation (|3.ip ] in that the 
the fluctuations become less, rather than more, anisotropic as k± increases. Isotropy is 
reached when k± ~ fcy ~ ki, where 

h ~ £ - 1 /2^3/2 _ fcoRo -2 (4 4) 

(k^ 1 ~ I Ft is sometimes called the Rossby deformation scale). At this wave number, 
the velocity is u(fcj) ~ (e/fl) 1 / 2 (using equations (|1.3p and l|1.4p ) and so the Rossby 
number at the corresponding scale is kiu(ki)/il ~ 1. Therefore, at k > ki, rotation 
is irrelevant and turbulence is of the familiar isotropic Kolmogorov type, with E(k) ~ 
£ 2/3£,-5/3_ Thjg jg 5 f course, the physically inevitable outcome because unlike in the 
case of magnetised turbulence, which can feel the mean magnetic field at any scale, the 
hydrodynamic turbulence cannot feel the mean rotation rate if the local Rossby numbers 
are large. It is reassuring that the critically balanced cascade predicted by the general 
argument proposed here has naturally led to this Kolmogorov state. 

The cascade path and the resulting spectrum are sketched in figure [T] We have illus- 
trated the case discussed above, where energy is injected isotropically and in the weak 
turbulence regime. More generally, we expect that energy injected at a given (k±, k^) will 
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injection 

— 1 ► 



Figure 1. A sketch of cascade path and spectra for the rotating turbulence: both the case of 
injection at k± — k« = fco (and Ro <C 1) and that at k« «C k±_ = &2D are shown (see sections 14, II 
and !4.3|l . In the case of polarisation alignment, scalings shown in the two upper panels should 
be modified as explained in section l4~5l 



travel towards the CB path [equation (|4.3p ] followed by isotropic Kolmogorov cascade 
(as, e.g., shown in figure [T] for the case of quasi-2D injection, discussed in section |4~3|) . 
Note that both the relationship (|4.3p between the parallel and perpendicular scales and 
the the wavenumber (|4.4|) of the transition to isotropy depend only on e and CI, but not 
on on the injection scale(s). Thus, the CB-to- isotropy path represents the "natural" state 
of rotating turbulence. We think that equation (|4.3[) should prove to be a good predic- 
tion for the relationship between the parallel and perpendicul ar correlation length s in 
the columnar vortical structures observed in experiments (e.g.. iDavidson et aLll2006h . 

Our predictions of the spectral slopes (|4.1|) and (|1.4|l and the scaling of the transi- 
tion wave number (|4.2[) with Ro are experimentally and numerically verifiable (but see 
section l4~5l for possible alternative scalings). A change of spectral slope from ~ —2.2 
(perhaps consistent with —5/2) to ~ —5/3 appears to have been observe d in rotating 
turbulence experiments with small initial Ro (Fig. 5b in Morize et al. 20051) , although it 
is premature to say if these results are definitely related to the transition to CB or merely 
to instrumental noise at high wave numbers. What does seem to be known definitely is 
that rotating turbulence has a clear tendency to a state where the local effective value 
of the Rossby numbe r is ~ 1, i.e., the linear and nonlinear time scales are comparable 
( Davidson et q/Jfeood ). This is consistent with the transition to CB that we are proposing 
and the nontrivial prediction is that while the spectrum is Kolmogorov for wave numbers 



Critical balance: towards a universal scaling conjecture 



7 



above this transition k± > k± c [equation (|1.4)l ]. the turbulence remains anisotropic up 
to the second transition wave number (14.41). 



4.2. Local scale- dependent anisotropy 

Quantitatively checking equation (|4.3[) is nontrivial. One possibility is to measure the 
energy spectrum with respect to parallel wave numbers: by definition, k^E(k^) ~ u\ ~ 
k±E(k±), so, using equations (| 1 .4|) and (|4.3[) . we find a distinctive scaling: 

E(k\\) ~ e 4/5 n- 2/5 k~ 7/5 , k < k\\ < h, k ±c <k±<h (4.5) 

(note that this is only valid in the CB regime). 

There is a similar result for Alfven-wave turbulence based on equation (|3.1[) : 

E{k\\) ~ ev^K 2 , k\\ > k , k± > k ±c (4.6) 

which appears to have recently be en corroborated by the solar wind measurements 
(jHorburv et ~al\ 120081 : IPode sta 2009). It is from the MHD experience that one learns 



about an important subtlety in the definition of fen or, more precisely, of the parallel 
scale 1 1| ~ fc^ 1 [in practice, scalings are usually extracted via structure functions rather 
than spectra: Su 2 ± (l) ~ l~ 1 E(l~ 1 )]. In order for the scale-dependent anisotropy to be- 
come apparent, 2|| had to be defined with respect to the "local mean field," i.e., the global 
mean magn etic field plus its perturbations at all scales larg er than the one we are in- 
terested in ( Cho fc Vishniac l200d: iMaron fc Goldreich|[200lT . Physically, this is because 



an Alfvenic perturbation can only "see" the local field, not the globally averaged one. 
Mathematically, measuring lu with respect to the global mean field would only capture 
the anisotropy of the largest-scale perturbation, while for all smaller-scale ones, such 
globally defined In "slips off" one field line to a neighboring one and effectively picks up 
perpendicular variation rather than the parallel one. 

Similarly, for rotating turbulence, we anticipate that some scheme might have to be 
devised to measure parallel correlations along the local mean vorticity direction rather 
than along the global rotation axis. Indeed, it is physically intuitive that the inertial 
waves would propagate along the total vorticity u) = 2ft + Suj, where 8u) = V x u (see 
further discussion in appendix [A"|) . In the CB regime, this deviation, while significant for 
measuring fen, is small: Suj/Q, ~ ku/VL ~ u/v ~ kn/kj_ "C 1; once isotropy is restored, 
Suj/fl ~ 1 (so inertial waves no longer have a definite direction of propagation). 

4.3. Is inverse cascade possible in rotating turbulence? 

Another interesting experimental possibility would be to stir the turbulence in a 2D 
way and find out whether it will develop an inverse cascade, bringing it first to the CB 
state and then to the isotropic Kolmogorov state. This possibility depends on the inverse 
cascade proceeding at a rate larger than the 2D structures radiating inertial waves, with 
energy thus directly transferred into the CB state (see the argument at the end of section 
[2] regarding the instability of a 2D state). 

A putative inverse cascade followed by a direct critically balanced cascade is sketched 
in figure[TJ Suppose the energy is injected at k± = &2D ^ &||- The inverse energy cascade, 
if it is sustainable, will give rise to the spectrum (|1.4| for kj_ < /c2D- This will extend, 
presumably at constant fcii, to k± ~ k± c , given by the first expression in equation (|4.2|) . 
At this point the turbulence is again 3D and the cascade should "turn around" and follow 
the CB path as before. Interestingly, the net perpendicular energy flux (integrated over 
fcii ) is zero for k± < although the spectrum is fc ± J ^ 3 extending to wave numbers both 
larger and smaller than with the infrared cutoff given by k± c . Since the velocity at 
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k± c is u ims ~ £ 2 / 5 (fe||f2) -1 / 5 , we have k± c ~ (k\\p,/u ms ) 1 /' 1 , where fey is the parallel wave 
number at which the energy was injected. Note that the cascade reversal is a nontrivial 
consequence of anisotropy; in isotropic turbulence, zero flux would imply thermodynamic 
energy equipartition, E(k) oc k 2 . 

4.4. Numerical evidence and finite-box effects 

There is a large body of literature on numerical simulations o f rotating turbulence . The 
two most recent and best resolved numerical studies are byE lnmni et all (|2009l ) and 



Thiele fc Mullen ((2009). We refer the reader to these papers for a comprehensive list of 



references to previous numerical work, which we will not replicate here. Let us discuss 
the results, which appear to be consistent in many independent investigations. 

Relating numerical evidence to scaling theories like the one proposed above is far from 
straightforward because simulations are typically done in periodic boxes and we have not 
discussed the effects of finite dimensio ns of the contain ing volume on wave turbulence. 
For MHD turbulence, it was shown by iNazarenkol ( 2007 ) that a finite box size along the 



mean magnetic field can lead to suppression of nonlinear interactions between modes with 
different fey . This results in a qualitatively different evolution of the 2D non-propagating 
fcii = mode and the wave modes with finite fey . Simil ar effects are probabl y operative 
in the numerical simulations of rotating turbulence (cf. ICambon et a7] l2004h . especially 
in re latively shallow boxeju Both papers cited above (as well as many previous ones, 
e.g., ISmith fc Walemj I1999T ) report significant accumulation of energy in the fey = 
modes, via a nonlocal inverse cascade the fey = modes also have a different spectrum 



than the finite-fey modes. Note that in the local inverse cascade scenario of sect ion I4~3| we 
envisioned energy injection at very small, but finite k\\ and did not consider the dynamics 
of the the exact fey = modes (whose existence is particular to numerical boxes). 

Another feature of the numerical simulations where the Rossby numbers associated 
with the forcing scale are low (the case we are considering in this paper), is what ap- 
pears to be a robust fcj 2 scaling of the energy spectrum (see papers cited above and 
references therein) . Is this a contradiction with the scaling predictions of section 14.11 ? A 
similar problem was encountered in interpreting numerical simulations of MHD turbu- 
lence, which consistently give E ~ k ± 3 ^ 2 (M aron fc Goldreichl 2001 : Mason et al. 2008 ). 



5/3 i | 

rather t han k , that followed from the scaling arguments of section [31 To explain these 
results, iBoldvre v (200G) proposed a modification of the standard CB argument based on 



the idea that nonlinearity is weakened in scale-dependent fashion if the fluctuating fields 
align in a certain way. It turns out that a similar modification can be constructed for 
rotating turbulence and yields a spectrum that agrees with numerical evidence. 



4.5. Polarisation alignment 

The estimate (|1.2j) for the nonlinear decorrelation time was correct subject to assuming 
implicitly that fluctuations are not polarised in any particular way that might weaken the 
nonlinearity, i.e., that the direction of u± decorrelates over the same scale as its ampli- 
tude. Let us suppose instead that a typical turbulent fluctuation is thrce-dimensionally 
anisotropic with characteristic wave numbers k x S> k y 3> k z = fey , where x is the direc- 
tion of maximum gradients remaining approximately the same throughout the fluctuation 
and z the direction of the propagation of the inertial waves. Then, since in a system with 
fey fcj_ the perpendicular velocity is individually incompressible, Vj_ • u±_ = 0, we have 

f Because of the anisotropy, even cubic boxes are effectively shallow. T his is also known in 
MHD turbulence, where simulations are routinely done in long boxes (e.g., [Maron fc Goldreichl 
120011 : iMason et aZ.ll2008ft . 
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u y (k x /k y )u x S> u x . Note that if we took k y = 0, we would simply have a monochromatic 
inertial wave, which, as it is easy to show, is an exact nonlinear solution (see appendix 
[A"| . However, if a wave packet of such waves is introduced, there would be nonlinear 
interaction and we are now attempting to determine how much the inertial-wave-like 
polarisation of fluctuations can be preserved in a strongly turbulent nonlinear state. 
To estimate the nonlinear decorrelation time, we now replace equation (|1.2p with 

t nl ~ k v u v ~ k±uj_(k±)6(k±), (4.7) 

where kj_ ~ k x , u± ~ u y and 6 ~ k y /k x ~ u x /u y <C 1 is the velocity angle responsible for 
weakening the nonlinearity (9 — would correspond to an inertial wave). To determine 
this angle, we need a physical hypothesis about the degree to which the inertial-wave 
polarisation of the velocity field is preserved across a typical turbulent fluctuation. We 
argued in section l4~2l that . physically speaking, inertial waves should propagate along the 
local mean vorticity direction rather than the global rotation axis and so the direction 
of anisotropy is scale-dependent. Thus, all directions within a fluctuation are determined 
to within an angular uncertainty 69 ~ 6uj±/£1 set by the local value of the perpendicular 
vorticity fluctuation. One might then postulate that that 9 ~ 89, i.e., that there is 
a tendency to preserve the inertial-wave polarisation to the maximal possible degree 
(polarisation alignment conjecture). Then 

e ~ Hull ~ h^± ^^jL^h. (4 8 ) 

Q Q v ky 1 

where we have taken u\\ ~ u y (as in an inertial wave), and used the CB conjecture: 
t nl ~ k|| u > where v ~ fl/k x and tnl is given by equation (14. 7|) . Since 8 ~ k x /k y by 
definition, we find k y ~ (k\\kx) 1 ^ 2 , 9 ~ (fcy/fc^) 1 / 2 , and so equation (|4.7p becomes 

tnl 1 ~ (*x*||) 1/2 «±(fe±). (4.9) 
Now combining this with equation (|1.3p and assuming CB, r ~ tnl ~ (k\\v) _1 , we find 

£7(&l) - [ev(fc ± )] 1/2 fc^ 3/2 - (eO) x / 2 A;- 2 , (4.10) 
fc|| ~ s 1 ' 2 [v(k ± )}- 3/2 k 1 / 2 ~ E Van-3/a fe a _ (4 u) 

These formulae replace equations (|1.4[) and (|4. 3|) for rotating turbulence with polarisa- 
tion alignment. Equation (|4.f f p implies that isotropy is again restored at the Rossby 
deformation wave number [equation (|4.4|) ]. while the transition wave number from weak 
to critically balanced turbulence is, instead of equation (|4.2p . determined by 

fc±c ~ £" 1/4 fej /2 3 / 4 ~ fcoRo" 1 . (4.12) 

The sketch of the cascade path in figure Q] is still valid, but with the new scalings for 
the CB regime (so at k± = k± c , the spectral slope now changes from —5/2 to —2 and 
at k± = ki, from —2 to —5/3; for k± c < k± < ki, ku cx k 2 ^). The parallel spectrum is 
derived as in section 14.21 and so equation (|4.5p is replaced by 

E(k\\) - e [vikxT 1 k^ 2 - e 3/4 fT 1/4 fc[ 3/2 . (4.13) 

The argumen t presented above is more or less analogous to the argument proposed by 



Boldvrev (2006) for MHD turbulence: he conjectured polarisation alignment between the 



perpendicular velocity and magnetic field fluctuations, which amounts to assuming that 
an Alfven-wave polarisation is approximately preserved. The scalings he derived can be 
read off from equations (|4.10p . (|4.1ip and (|4.13p by setting v(k±) = va instead of 2£!/fcj_ 
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(note that equation (14. 6|) remains unchanged). Numerical simulations co nfirmed these 
scalings, including also th e scale-dependent alignment between the fields ([Mason et all 
2008; Bo ldvrev et af]|2009h . It appe ars plausible that the k , 2 spectra observed in simula- 
tions of forced rotating turbulence ([Mininni et al. 1 120091: iThiele fe Muller|[2009h could be 
similarly explained by the polarisation alignment argument we have given here. However, 
a word of caution: most solar wind measurements of Alfvcnic turbulence do not support 

3/2 5/3 i I I 

the k~ scaling and incline rather towards k]_ (jBale et all 2005; H orburv et aLll2008 ; 
Sahraoui et al. 2009). Thus, it remains unclear whether polarisation alignment occurs in 



nature as well as in numerical boxes. We stress that the possibility of polarisation align- 
ment does not undermine CB as a universal scaling conjecture — what is revised in this 
version of CB turbulence is the assumption of isotropy in the perpendicular plane. 



5. Stratified turbulence 



Another hydrodynamic example where a critically balanced cascade should emerge is 
the stably vertically stratified turbulence, anisotropic with fc || ^> k \ , whe re k\\ and k± are 
the vertical and horizontal wave numbers, respectively ( Laval et aZ.ll2003f ). The dispersion 



relation for (incompressible) gravity waves is uj = Nk±/k, where N is the Brunt -Vaisala 
frequency. Since the roles of fc|| and k± are reversed compared to the MHD and rotating 
turbulence, the arguments presented above have to be somewhat modified. 

First, the incomprcssibility V • u = now implies U|| ~ (k±/k\\)u± <C uj_ (since 
fcj_ <C fcii, it is no longer true that Vj_ • u± = as was the case for MHD and rotat- 
ing turbulence). This implies that the nonlinear interaction time continues to be given 
by equation (|1.2p . If CB is assumed, the horizontal energy spectrum is, therefore, still 
Kolmogorov [equation (|1.4|) ]. while the relationship between the horizontal and vertical 
wave numbers is [analogous to equations p.ip and (|4.3p ] 

k± ~ eA~ 3 fcf = l 2 k\, (5.1) 



where lo = e l / 2 N 3 / 2 is called the Ozmidov ( 19921 ) scale. Using equation ()5.1|) . we can 



calculate the vertical energy spectrum corresponding to the horizontal spectrum (| 1.4)1 in 
a way analogous to the derivation of equation 



E{k\\) ~ N k7 (5.2) 



(this spectrum was prey iously proposed on dimensional grounds by iDewanl 1997 and by 



Billant fc Chomazll200lh . This argu ment is basically a reformulation in the CB language 



of the scaling hypothesis proposed bv lLindbord ( 2006t h where we also refer the reader for 



discussion and references on the atmospheric measurem ents and numerical simulation s, 



which appear to be consistent with this theory (see also iLindborg fc Brethouwerll2007h . 

The situation here is similar to rotating turbulence: the anisotropy gets weaker as fc|| 
increases. The cascade becomes isotropic at the Ozmidov scale, kj_ ~ fen ~ Iq , where 
the local Froude number u±/loN ~ 1. At smaller scales, the turbulence cannot feel 
the mean gradient and becomes isotropic. Then both horizontal and vertical spectra are 
Kolmogorov, so there should be a spectral break at the Ozmidov scale for the vertical 
spectrum [transition from (|5.2[) ] . but not the horizontal one, which is Kolmogorov already 
in the CB regime. The cascade path is similar to figure [1] with fc|| and k± swapped and 
equation (|5.ip used for the CB line. 
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6. Discussion 

We have proposed that the critical bala nce of linear and nonlinear time scales, orig- 
inally introduced for Alfvenic turbulence (|Higdonlll984l : iGoldreich fe Sridharlll995[) and 
more recently, for other types of magnetised plasma turbulence ( Cho fc Lazarianl 12004 
ISchekochihin et al. 20091 ). should be used as a universal scaling conjecture for anisotropic 
wave systems. While t here are some indications that this idea works for stratified turbu- 
lence |Lindbord 1200(1 . it has not been tested in rotating turbulence, for which we have 
proposed a novel energy cascade scenario and a set of verifiable predictions. 

In neutral fluids, the two examples we have considered — rotating and stratified tur- 
bulence — suggest that the critically balanced cascade provides a path from the strongly 
anisotropic fluctuations caused by the presence of the external field (mean rotation or 
gradient) to isotropic Kolmogorov turbulence inevitable at sufficiently small scales (cf. 



Davidson et al. 20061) . In that way, neutral fluids are different from conducting fluids 



and plasmas, where the presence of a mean magnetic field is felt at all scales and the 
anisotropy only gets stronger as smaller scales. 

We observe that one might draw parallels betw een the CB principle and the "gener- 
alised Phillips spectra" ( Newell fc Zakharov 20081 ) that are thought to emerge from the 
break down of weak turbulence in ma ny wave system s, e.g., surfac e water waves (IPhillipsl 
1958lll985tlNewell fc Zakharovteooll . Rossby waves (|Rhinedll975h. Kelvin waves in crvo- 



genic turbulence ( Vinenl 20001 ) and Bose-Einstein condensates dProment et al. 20091). See 
also discussions of such universal flux-independent spectra by Denissenko et all (|2007 ) 
and ICardv et al. (I2008L p. 12). Some of these examples do not quite fit the anisotropic 
type considered here and the break down of weak turbulence there may be related to 
formation of singular structures: this is discussed in more detail in appendix [Bj 

We acknowledge important discussions with P. Bcrloff, S. Cowley, P. Davidson, J. 
McWilliams, F. Moisy and A. Newell. We are also grateful to the anonymous referees for 
raising a number of issues, which led to significant expansion and improvement of the 
paper. This work was supported by STFC (A.A.S.) and the Leverhulme Trust Network 
for Magnetised Plasma Turbulence. 



Appendix A. Reduced equations 

A remarkable simplification of the underlying dynamical equations for turbulent fluc- 
tuations can be achieved by systematically taking into account their anisotropy. Here 
again CB serves as a guiding principle, but this time as an ordering assumption: the 
expectation that the linear and nonlinear time scales would be comparable tells one how 
to order the fluctuation amplitudes with the expansion parameter e = fc||/fcj_. This leads 
to reduced systems of equations, which are often more transparent physically and require 
less computational power to simulate numerically. An additional advantage of simulating 
reduced equations is that the transition to the asymptotic anisotropic regime is carried 
out analytically and so does not eat up resolution. 

Two well known examples of such reduced systems are the Reduced MHD (RMHD) 
equations for the Alfvenic turbulence (reproduced at the end of this section) and Elec- 
tron Reduced MHD (ERMHD) equations for the turbulence o f kinetic Alfven waves at 
sub-Larmor scales (both are derived under the CB ordering in Schekochihin et a/1l2009[ 
but the RMHD has been known since Strauss! 19761 while equations mathematically sim- 
ilar to ERMHD hav e been used to describe whistler turbulence in plasmas; see, e.g., 
Biskamp et aZ.I fl999h . A kinetic reduced system that em erges from the same principles 



is gyrokinetics (jFrieman fc Chenlll982t iHowes et aLll2006l ). a general description of mag- 
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netised plasma turbulence of which RMHD and ERMHD are particular limits. All of 
these equations have been successfully s imulated numerically (jPerez fc Boldvrevl 12008 : 
Cho &: Lazariaij|2009t iHowes et al\\200$) and in the case of RMHD and ERMHD, the 



results are known to be asymptotically consistent with simulations of the unreduced 
equations. 

Here we explain the procedure for deriving a reduced system on the example of rotating 
turbulence, showing again its very close resemblance to magnetised fluid systems and 
providing concrete justification for some of the assumptions made in the main text. 

Our starting point is the Navier-Stokes equation for an incompressible fluid of viscosity 
v and density p = 1 rotating at the rate CI = Qz (z is the rotation axis): 

du 

— +u- Vit + 20 x u = -Vp + v\J 2 u, (Al) 

where p is pressure determined by V • u = 0. This equation supports inertial waves with 
frequency u = ±2Qk^/k and corresponding eigenfunctions u — (±ik z /k,l,^fik x /k)u v , 
where k — (k x ,0,k z ) without loss of generality. Note that for an inertial wave, the 
perturbed vorticity is aligned with velocity, 5u> = =pku, and so a monochromatic inertial 
wave is an exact nonlinear solution of equation (|A 1[) (because u ■ Vit = Su) x u + Vit 2 /2; 
the gradient part can be absorbed into pressure). 

In the anisotropic regime (low Rossby numbers), equation (|A ip is expanded in a small 
parameter e = k\\/kj_- Using the CB as an ordering prescription, we order the time 
scale of the fluctuations as u> ~ ed ~ u± ■ V j_. We also order u\\ ~ u± guided by the 
relationship between them in an inertial wave. 

To lowest order in e, V-u = - u±_ = 0, so the perpendicular motions are individually 
incompressible and can be represented by a stream function: u±_ = z x Vx$. In the next 
order, the incompressibility condition allows us to find the divergence of the second-order 
correction to u± (to be useful shortly): 

V • u = V± ■ m ( , 2) + ^ = 0. (A 2) 

oz 

The perpendicular part of equation (|A 1[) is (keeping two lowest orders) 

+ u± ■ V±u± - v^7\u± = -2tiz xu±- V±p. (A 3) 

In the lowest order, the left-hand side disappears and so right-hand side must vanish too. 
This gives p = 20$. Now taking the perpendicular curl (Vix) of equation (|A3|) . we get 

|-Vi$+{$,Vl$} =2Q^+i/Vi$, (A 4) 

at L ' oz 

where {f,g} — z ■ (Vj_/ x Vj_^) and we have used equation (|A 2|) to express x 

(2flz x ) — 2Q.Z, Vj_ • u^ 1 . Finally, taking the parallel part of equation (|A ip to lowest 

order and using p = 2f2$, we get 



^ + {*,«„} =-20 1| +I/VI*,,. ( A 5 ) 

Equations (|A 4p and (|A h\ are the desired reduced system, which we will refer to 
as Reduced Rotating Hydrodyna mics (RRHD) . Up t o notational differences, they are 
equivalent to the reduced system of Julien et al. ( 19981 ) (derived under slightly differently 
formulated assumptions). Let us itemise some of their properties: 

(a) they support inertial waves with lo = ±2ttk^/k± and corresponding eigenfunctions 
Uii = ±fc^$ (so, the wave is circularly polarised: as it propagates along z, the velocity 



Critical balance: towards a universal scaling conjecture 13 

vector rotates in the plane perpendicular to k±); inertial-wave packets with fixed k± and 
arbitrary amplitude are exact nonlinear solutions; 

(b) the velocity is u = z x Vi$ + zuu and the (perturbed) vorticity 5us = — z x 
V_lU|| + zV\^\ therefore, the RRHD equations can be written (omitting viscosity) 

— Vi$ + {$,Vi$} = (5U + 2SI)- V«||, = -{5uj + 2Q.) • V$, (A 6) 

thus, besides the 2D self-advection of the perpendicular velocity, the basic physical pro- 
cess is propagation of inertial waves along the total vorticity lines (recall section [472]) ; 

(c) the 3D nature of the turbulence is enforced via linear propagation terms (recall 
section [2]); when they are present, the system conserves one invariant, kinetic energy 
/ d 3 r (| Vj_<fr| 2 /2 and should have a direct cascade; when d/dz — 0, there are three 
invariants: perpendicular kinetic energy J d 3 r\ V^$| 2 /2 (inverse cascade; recall section 
14. 3p . enstrophy J d 3 r | V 2 L $| 2 /2 (direct cascade), and parallel kinetic energy J d 3 ru 2 /2 
(direct cascade of a passive quantity) ; 

(d) while the RRHD equations were derived under the CB ordering, they allow both 
the weak and the strong turbulence regimes and so can track the transition from the 
former to the latter (section [4. 1[) ; they also remain valid if polarisation alignment occurs 
(section [4.51) and so can be used to measure and study it. 

(e) when nondimcnsionalising the RRHD equations, one can rescale the parallel and 
perpendicular distances separately (subject to appropriate rescaling of the amplitudes), 
i.e., the aspect ratio of the "box" is formally infinite — this is because the anisotropy 
of the fluctuations was the basis for the asymptotic expansion that led to RRHD; note 
that the scaling arguments of section 0] imply that the anisotropy of rotating turbulence 
diminishes with scale, so equations (IA 4p and (|A 5[) will produce solutions that, at suffi- 
ciently small scales, violate the ordering assumptions behind the equations — this should 
be manifested by the development of ever finer parallel structure (see section 14. 2|) . 

Finally, for comparison, let us give the RMHD equations for the Alfvenic turbulence 
(derived from MHD in ex actly the same way as RRHD was from equation (|A 1[) : see, e.g., 
Schekochihin et al. 2009, section 2) 



^vi$ + {$,vi$} =w A ^vi* + {*,vi*} =v A h- vvi*, (A7) 

<9<J d$ 

-t-=v a — + {■$,$} = v A b-V$, (A 8) 

at oz 

where $ and VP are the stream functions of the perpendicular velocity and magnetic 
fluctuations: itj_ — z x Vj_$, 5B± = z x Vj^ (unlike in the rotating case, the parallel 
velocity and magnetic field fluctuations decoup le and are passive with respect to the 



perpendicular ones; see Schekochihin et al. 2009), and b — z + SB±/Bq is the direction 



of the total magnetic field (along which the Alfven waves manifestly propagate; recall 
section [4.2[) . The similarities with the equations for rotating turbulence are obvious, but 
there are also differences originating from the physical differences between the inertial 
and Alfven waves: the latter are non-dispersive (u = ±kuv A ), the perpendicular velocity 
and magnetic field fluctuations are linearly polarised and the eigenfunctions <f> = ±^ 
represent exact nonlinear solutions for arbitrary-shaped wave packets. 



Appendix B. Critical balance and coherent structures 

It might appear tempting to relate the CB principle to the emergence of coherent 
structures for the following two reasons. Firstly, some coherent structures arise due to 
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wave breaking, like in the system of water surface gravity waves or internal gravity waves, 
and such a wave breaking occurs precisely when the nonlinearity becomes of order of the 
linear contributions. Secondly, in some well known examples of coherent structures, like 
solitons or shocks, the linear and the nonlinear terms are in balance. 

First, consider the coherent structures that result from wave breaking. Such structures 
are typically singular, e.g., the wave profile has a derivative discontinuity. To be specific, 
consider the water-surface gravity waves, where the CB approach, i.e. , scale- by-scale 
balance of the linear and the nonlinear time scales, gives the well-known IPhillipsI (|l958h 
spectrum. Its connection to sin gular wave breakin g structures has been widely discussed 
since it was first suggested by Kadomtsev ( 1965 ) and later adopted by IPhillipsI ( 19851 ) 
himself (his original 1958 paper does not mention of wave breaking). However, there 
is an unce r tainty related to the geometry and the dimensionality of the wave crests. 
Kuznetsovl ( 2004 ) argued that the Phillips spectrum should correspond to wave crests 



with singularity at isolated points (cone-like shapes) whereas more realistic ID crests 
would give a different spectral exponent. Furthermore, coherent structures of different 
strengths or sizes can coexist and the resulting s pectrum can d epend on their distribution 
(e.g., in his refinement of the original theory Phillips! 1 19851 introduced a distribution 
function for the crest lengths per unit area of the water surface). Thus, there is no obvious 
universal link between the CB state and singular coherent structures of the wave-breaking 
type: there may be structures with spectra different from the CB prediction, but one can 
also imagine a CB state without any singular structures at all. 

Now let us turn to the nonsingular coherent structures. The most basic of the relevant 
nonlinear models is the Korteweg-de Vries (KdV) equation, u t + uu x + fiu xxx = 0, 
where /x is the dispersion coefficient. Does CB work for KdV turbulence? Naively, the 
idea might seem promising because, the KdV model predicts formation of solitons - 
coherent structures in which the nonlinear and the linear terms are balanced, in the 
spirit of CB. Balancing the nonlinearity and dispersion scale by scale (the second and 
the third terms in the KdV equation), we get E(k) ~ fc 3 . However, if the separations 
between the solitons are much greater than their widths, then for the scales intermediate 
between the soliton width and the inter-soliton separation, the spectrum is that of a 
set of delta-functions, so E(k) ~ const, which is very different from the CB prediction. 
This is because the interaction is nonlocal in scale, whereas CB assumes locality. Note 
also, that this is a ID dispersive system and the physical arguments in favour of CB in 
anisotropic wave-supporting environments presented in section do not generalise to it. 

In conclusion, there does not seem to be a universal link between the CB and coherent 
structures. In some systems there may be coherent structures but not CB because the 
latter assumes locality, which is not automatically guaranteed. Even if both coherent 
structures and a CB spectrum are present, the former need not be the cause of the 
latter as coherent structures might occupy a negligible volume. Finally, we reiterate that 
the physical argument for the formation of a CB state in MHD, which, as we showed 
above, may similarly be applied to the rotating turbulence, does not invoke either wave 
breaking or coherent structures and is rather based on local nonlinear energy transfer 
and anisotropic spatial decorrelation arguments (section [2]) . 
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